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PRODUCING  DESIRED  COLOR  STIMULI  WITH 
TRICHROMATORS:  A  TUTORIAL  ON  TECHNIQUE 


1.  SUMMARY 

This  report  summarizes  the  theory  of  color  mixing  when  using  a 
trichromator  system.  The  four-channel  colorimeter  system  described  is  of 
a  generic  nature.  The  derivation  of  formulae  needed  to  compute 
chromaticity  coordinates  and  luminance  requirements  are  presented.  A 
numerical  example  is  provided  which  demonstrates  the  application  of  the 
formulae. 

2.  INTRODUCTION 


The  need  sometimes  arises  in  laboratory  experimentation  to  produce 
stimuli  of  precise  chromaticity  and  luminance.  The  need  this  author  has 
encountered  involves  color  matching  experiments,  where  two  stimuli  are 
generated  that  must  be  matched  along  one  dimension  (either  dominant 
wavelength,  excitation  purity,  or  luminance).  In  the  author's  case,  the 
two  stimuli  formed  a  bipartite  field  (Greene  &  Laycock,  1984).  One-half 
of  the  matching  field  was  set  by  the  experimenter  to  be  a  predetermined 
and  fixed  chromaticity  and  luminance.  The  subject  adjusted  the  other 
half  of  the  bipartite  field  for  a  match  in  one  of  the  above  mentioned 
dimensions. 

To  produce  the  two  stimuli  for  color-matching,  a  trichromator  system 
is  described  which  mixes  the  appropriate  wavelengths  in  proper  luminance 
proportions  to  generate  the  bipartite  field.  Although  this  report 
describes  a  generic  colorimeter,  it  is  representative  of  ones  currently 
existing  in  laboratories  in  the  Crew  Systems  Effectiveness  Branch,  Air 
Force  Aerospace  Medical  Research  Laboratory  (AFAMRL/HEF),  Wright- 
Patterson  Air  Force  Base,  Ohio,  and  in  the  Human  Factors  Group,  Royal 
Aircraft  Establishment,  Farnborough,  England. 

3.  TRICHROMATOR  DESCRIPTION 


To  generate  stimuli  for  color-matching  experiments,  a  trichromator 
system  was  chosen  by  this  author  to  produce  light  sources  of  different 
wavelengths  for  mixing.  Depending  upon  how  much  luminance  is  required  of 
the  stimuli,  the  light  source  for  the  trichromator  system  can  either  be  a 
single  source  divided  among  the  three  monochromators,  or  multiple  light 
sources,  one  for  each  of  the  three  monochromators. 

Because  of  the  exact  nature  of  the  mathematical  equations  used  in  the 
calculation  of  wavelengths,  it  is  desirable  to  have  monochromators  whose 
selectable  wavelengths  are  accurate  to  0.1  nm.  Otherwise,  the  mixing  of 
wavelengths  will  result  in  less  than  desired  precision  in  resultant 
chromaticity  coordinates. 
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With  a  trichromator,  monochromatic  light  from  the  three 
monochromators  is  combined  through  four  optical  channels  in  such  a  way  as 
to  produce  two  stimuli.  In  order  to  do  this,  one  of  the  monochromators 
is  fixed  at  an  anchor  wavelength,  usually  in  the  longer  wavelengths  of 
the  visible  spectrum  (i.e.,  600+  nm).  The  light  from  the  anchor 
wavelength  is  split  in  half  by  use  of  a  beam-splitter,  forming  two 
channels.  Light  from  one  of  the  other  two  monochromators  is  combined 
with  one  channel  of  the  anchor  wavelength,  in  exact  quantities  of 
luminance  from  each,  to  produce  one  stimulus  of  a  precise  chromaticity. 

In  the  same  way,  the  third  monochromator  combines  monochromatic  light  of 
a  specific  luminance  with  the  other  channel  of  the  anchor  wavelength  to 
produce  a  second  stimulus,  again  of  a  required  chromaticity. 

Figure  1  shows  diagrammatical ly  the  trichromator  system  presently 
configured  at  the  Royal  Aircraft  Establishment.  The  four  optical 
channels  are  labelled  for  ease  of  understanding.  In  this  example,  light 
from  channels  2  and  4  form  the  anchor  wavelength  channels.  The 
monochromator  in  channels  2  and  4  is  always  set  at  the  anchor  wavelength. 
The  light  output  from  channels  1  and  2  is  mixed  to  form  the  chromaticity, 
of  a  controlled  luminance,  required  for  one  of  the  stimuli.  Similarly, 
the  output  from  channels  3  and  4  is  combined  to  produce  the  chromaticity 
of  the  second  stimulus,  again  at  a  precise  luminance  level.  This 
combination  of  light  outputs  is  shown  in  Figure  1.  The  AFAMRL 
trichromator  system  is  nearly  identical. 

4.  COLOR  MIXING 


The  tristimulus  values  resulting  from  the  mixture  of  two  or  more 
lights  are  the  sums  of  the  lights'  individual  tristimulus  values.  For  a 
mixture  consisting  of  only  two  lights,  the  point  representing  the 
chromaticity  of  the  mixture  is  located  on  the  straight  line  connecting 
the  points  representing  the  chromaticities  of  the  two  lights  on  one  of 
the  Commission  Internationale  de  I'Eclairage  (CIE)  chromaticity  diagrams 
(Graham,  1965,  pp.  270-294).  Any  of  the  CIE  chromaticity  diagrams  may  be 
used  for  the  purposes  described  herein,  but  the  present  discussion  will 
deal  only  with  the  1931  CIE  diagram,  as  it  is  the  most  familiar.  It  is 
this  fact  of  color  mixing  which  describes  the  basis  for  the  theory  of  the 
trichromator  system.  On  the  straight  line  drawn  between  two  points  on 
the  spectrum  locus  of  the  1931  CIE  diagram  lie  all  the  chromaticity 
coordinates  achievable  by  the  mixture  of  these  two  lights. 

There  are  three  steps  this  author  went  through  to  use  the 
trichromator  system  for  color  mixing.  The  first  was  to  select  the  anchor 
wavelength  for  mixing.  It  was  the  author's  desire  to  produce  stimuli  at 
a  fairly  high  luminance  level  -  2000  cd/m2.  It  was  this  consideration, 
coupled  with  the  desire  to  examine  dominant  wavelengths  in  the  longer  end 
of  the  visible  spectrum  (640  nm),  which  influenced  the  choice  of  the 
anchor  wavelength.  The  luminance  calibration  of  the  light  output  from 
the  monochromator  in  channels  2/4  revealed  the  need  for  an  anchor 
wavelength  of  650  nm  for  the  particular  color-matching  experiment  of 
interest. 
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Figure  1.  Optical  Representation  of  the  Trichrimator  Hardware. 
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The  second  step  was  to  calculate  the  variable  wavelength  for  mixture 
with  the  anchor  wavelength.  This  variable  wavelength  determination  is 
dependent  upon  what  stimulus  chromaticity  required. 

The  third  step  was  to  calculate  the  luminance  proportions  of  the 
anchor  and  variable  wavelengths  required  to  produce  a  given  stimulus 
chromaticity  at  a  precise  luminance.  The  next  two  sections  describe  the 
calculations  required  to  achieve  steps  2  and  3.  The  anchor  wavelength 
setting  can  vary  depending  upon  the  experimental  conditions  under 
investigation. 

4.1  Computation  of  Variable  Wavelength 

Figure  2  illustrates  the  mixing  of  wavelengths  within  a  chromaticity 
diagram.  Point  "A"  on  the  diagram  shows  the  location  on  the  spectrum 
locus  of  the  chosen  anchor  wavelength*  650  nm  in  this  example.  Point  "B" 
represents  the  desired  chromaticity  within  the  diagram.  In  order  to 
determine  which  wavelength  must  be  combined  with  the  anchor  wavelength  to 
produce  the  desired  chromaticity,  the  following  procedure  is  used.  A 
straight  line  is  drawn  from  the  anchor  wavelength  on  the  spectrum  locus 
(Point  "A")  through  the  desired  point  within  the  chromaticity  diagram 
(Point  "B"),  and  extended  to  where  the  line  intersects  the  spectrum 
locus,  493  nm  in  this  example  (Point  "C")*  The  dominant  wavelength 
represented  by  Point  "C"  is  what  is  required  for  mixing  with  the  anchor 
wavelength  to  produce  the  desired  chromaticity. 

It  is  possible  to  produce  any  desired  chromaticity  within  the  diagram 
by  the  combination  of  twc  dominant  wavelength  sources  with  appropriate 
luminance  contributions  from  each  source.  The  mixture  of  the  two 
wavelengths  in  appropriate  quantities  determines  what  stimulus 
chromaticity  is  produced.  In  this  case,  the  appropriate  quantities  of 
monochromatic  light  to  be  mixed  are  expressed  in  terms  of  luminance  of 
each  wavelength.  The  calculation  of  the  amounts  of  luminance  required  at 
each  wavelength  is  discussed  in  a  later  section. 

The  graphical  method  described  above  is  not  only  tedious,  but  also 
does  not  yield  exact  answers.  To  circumvent  this,  a  mathematical 
solution  may  be  implemented  on  a  digital  computer. 

Trigonometry  is  used  to  compute  a  table  of  angles,  theta,  formed  when 
drawing  a  line  from  the  anchor  wavelength  to  every  point  on  the  spectrum 
locus  at  one-nanometer  intervals.  The  chromaticity  coordinates  of  points 
on  the  spectrum  locus  are  provided  at  one-nanometer  intervals  in  texts 
such  as  Wyszecki  and  Stiles,  1982,  pp.  725-735.  Figure  3  illustrates 
this  trigonometric  solution.  Figure  4  shows  graphically  the 
determination  of  theta  (referenced  to  an  anchor  wavelength  of  650  nm)  for 
a  specific  example  -  490  nm. 

A  copy  of  a  FORTRAN  program  written  to  determine  theta  angles  is 
found  in  Appendix  A.  From  this  FORTRAN  program,  once  the  anchor 
wavelength  has  been  selected  and  its  chromaticity  coordinates  entered,  a 
table  of  thetas  is  compiled  for  each  one-nanometer  point  around  the 
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Figure  2.  Trichromator  Operation  in  1931  CiE  Chromaticity  Space. 
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Figure  3.  New  Axis  Coordinate  Grid  for  Determination  of  Theta  Angie. 
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spectrum  locus  (360-830  nm),  referenced  to  the  anchor  wavelength.  A 
compilation  of  thetas  for  an  anchor  wavelength  of  650  nm  is  shown  in 
Appendix  B. 

For  the  sake  of  convenience,  this  tabulated  set  of  thetas  is  stored 
on  a  computer  as  a  data  file.  A  second  computer  program  then  can  be 
written  to  search  through  the  table  of  thetas  to  find  the  closest  integer 
wavelengths,  and  then  perform  the  necessary  linear  interpolation  to 
arrive  at  the  exact  wavelength  needed  for  mixing.  A  copy  of  a  BASIC 
program  written  to  achieve  this  purpose  is  found  in  Appendix  C. 

This  section  has  summarized  a  mathematical  solution  to  the  problem  of 
determining  the  wavelength  required  for  mixing  with  an  anchor  wavelength 
to  achieve  a  desired  stimulus  chromaticity.  The  next  section  describes  a 
procedure  for  determining  the  luminances  which  are  required  to  produce 
the  desired  chromaticity  at  a  specific  luminance. 

4.2  Luminance  Calculations 

With  the  trichromator  system,  the  luminance  contributions  from  the 
three  monochromators  is  controlled  by  a  series  of  calibrated  neutral 
density  wedges.  The  author's  practice  has  been  to  store  this  luminance 
calibration  data  on  a  computer  as  a  data  file. 

The  typical  situation  the  researcher  generally  faces  is  one  of  trying 
to  produce  a  specified  chromaticity  at  a  controlled  luminance  level. 

Once  the  anchor  wavelength  is  chosen,  two  more  steps  remain.  The  first 
is  to  determine  the  variable  dominant  wavelength  to  produce  the  desired 
point  on  the  chromaticity  diagram.  This  step  was  outlined  above  in 
Section  4.1.  The  second  step  is  to  calculate  the  luminances  required  for 
the  two  dominant  wavelengths  to  arrive  at  the  needed  luminance  level. 

In  order  to  accomplish  these  two  steps,  the  following  equations  are 
provided  and  were  derived  from  Hunt  (1975). 

Given  that: 


(L2  +  L2) 


and  that  LI  +  L2  =  L3  (3) 
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where: 


(XI, Yl)  =  CIE  chromaticity  coordinates  for  monochromator 
1  wavelength  setting  (channel  1), 

(X2,Y2)  =  CIE  chromaticity  coordinates  for  the  anchor 
wavelength  setting  (channel  2), 

(X3,Y3)  =  CIE  chromaticity  coordinates  for  the  desired 
point  inside  chromaticity  space, 

Ll  =  luminance  from  channel  1, 

L2  =  luminance  from  channel  2, 


and  L3  =  luminance  of  point  of  desired  chromaticity. 


Algebraic  manipulation  of  the  above  equations  results  in  the 
equations  for  the  luminance  necessary  from  each  monochromator.  The 
results  of  the  algebraic  manipulations  are  as  follows: 


(4) 

(5) 

where: 

c  =  ((Y3-Y1)  /  Yl) 
and  d  =  ((Y3-Y2)  /  Y2) 

5.  EXAMPLE 


Ll  =  {(-d  /  (c-d))  X  L3 
and  L2  =  ((c  /  (c-d))  X  L3) 


It  may  be  helpful  to  see  a  computational  example  using  the  above 
derived  formulae.  For  this  example,  it  is  desired  to  produce  a  stimulus 
whose  chromaticity  coordinates  are  x  =  0.2,  y  =  0.7  at  a  luminance  of 
500  cd/m2. 

The  anchor  wavelength  for  this  example  is  650  nm  (x  =  0.72599, 
y  =  0.27401).  The  first  task  is  to  draw  a  straight  line  from  the  650  nm 
point  on  the  spectrum  locus  of  the  1931  CIE  diagram,  through  the  desired 
chromaticity  (x  =  0.2,  y  =  0.7).  Figure  5  illustrates  this  step. 
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Figure  5.  Finding  Theta  Angie  and  Spectrum  Locus  intersection 
for  Numerical  Example. 
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The  next  task  is  to  calculate  the  resulting  theta  angle.  To  avoid 
the  ruler-and-protractor  approach,  the  aforementioned  computer  programs 
were  used  to  determine  the  theta  angle  of  140.997  degrees.  Examining  the 
data  table  of  theta  angles  provided  in  Appendix  B  reveals  the  required 
monochromator  setting  needed  for  mixture  with  650  nm  is  between  516  and 
517  nanometers  (between  141.3  and  140.7  degrees).  A  simple  linear 
interpolation  is  performed  to  find  the  monochromator  setting  needed  for 
mixture  with  the  anchor  wavelength.  This  calculation  produces  a 
wavelength  of  516.3  nm. 

Equations  4  and  5  define  the  luminance  requirements  from  each  of  the 
two  dominant  wavelengths  in  order  to  produce  the  required  stimulus  at  a 
given  luminance  level.  In  this  example,  the  desired  chromaticity  of 
X  =  0.2,  y  =  0.7  is  to  have  a  luminance  of  500  cd/m2. 

Upon  examining  these  equations,  it  can  be  seen  that  the  variables 
needed  for  input  are  the  (x,y)  chromaticity  coordinates  for  the 
following:  anchor  wavelength,  variable  wavelength  setting  for  mixture 
with  the  anchor  wavelength,  and  desired  chromaticity  within  the  1931  CIE 
diagram.  At  this  point,  two  of  the  three  pairs  of  chromaticity 
coordinates  are  known  -  the  anchor  wavelength  and  the  chromaticity  of  the 
desired  point  within  the  chromaticity  diagram. 

A  linear  interpolation  must  be  performed  on  the  x  and  y  coordinates 
(between  the  chromaticity  coordinates  for  516  and  517  nm)  to  obtain  the 
final  pair  of  coordinates  for  the  equations.  Linear  interpolation 
produces  the  resultant  coordinates:  x  =  0.04739,  y  =  0.82112. 

To  use  Equations  4  and  5,  we  define  our  variables  as  follows:  let 
XI,  Yl  be  the  coordinates  of  the  variable  monochromator  in  channel  1,  XI 
=  0.04739,  Yl  =  0.82112,  (516.3  nm  in  this  example).  Let  X2,  Y2  be  the 
coordinates  of  the  anchor  wavelength,  X2  =  0.72599,  Y2  =  0.27401,  (650  nm 
in  this  example).  Let  X3,  Y3  be  the  coordinates  of  the  desired 
chromaticity,  X3  =  0.2,  Y3  =  0.7. 

Substituting  values  into  Equation  4  we  have: 


-(0.7  -  0.27401) 

■  (Tl&TTOT - 

LI  =  _  X  500 

0.7  -  0.82112  -  0.7  -  0.27401 

- 0:82112  0.27401 


(6) 


Solving  Equation  6,  we  find  that  LI,  the  luminance  contribution  at  516.3 
nm,  is  456.8  cd/m2.  ^Thus,  via  Equation  3,  the  luminance  required  at  650 
nm  (L2)  is  43.2  cd/m^. 
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6.  CONCLUSION 


The  numerical  example  worked  in  Section  5  is  for  producing  one 
stimulus  of  a  required  chromaticity  at  a  desired  luminance  level.  The 
identical  procedure  would  be  followed  to  determine  the  other  dominant 
wavelength  needed  to  mix  with  the  anchor  wavelength  to  produce  a  second 
stimulus.  In  the  same  fashion  as  above,  the  luminance  contributions 
would  also  be  determined. 

The  theory  of  color  mixing  using  a  tri-stimulus  colorimeter  to 
produce  two  stimuli  of  known  chromaticity  and  desired  luminance  has  been 
presented.  The  derivation  of  the  required  formulae  involve  a  simple 
algebraic  manipulation  of  a  set  of  equations.  Since  this  manipulation 
has  not  heretofore  been  published  in  the  literature,  it  is  hoped  that 
the  sharing  of  this  basic  information  will  prove  useful. 
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Appendix  A  -  FORTRAN  Program  to  Compute  Theta  Angles 
FILENAME  ANCHOR. FTN 

GIVEN  AN  "ANCHOR"  POINT.  DRAW  A  CHORD  TO  EACH  OF  471  POINTS  ON 
1931  CIE  CHROMATICITY  DIAGRAM.  MEASURING  ANGLES  ACC0RDTN6  TO 
THE  FOLLOWING  CONVENTION: 

90 

180  ++  0 
270 


DIMENSION  X(471).  Y(471).  AN6LE(471) 

DATA  0/57.2958/ 

INPUT  471  PAIRS  OF  (X.Y)  COORDINATES  DEFINING  CHROMATICITY 
LOCUS  FROM  DISK  FILE 

OPEN  ( UNIT=1 . NAME= ' CMF31 . DAT ' . TYPE= ' OLD ' . ACCESS^ ’ SEQUENTIAL ' ) 
READd.*)  (Xd).  Yd).  I  =  1.  471) 

CLOSE  (UNIT=1) 

INPUT  X.Y  COORDINATES  OF  ANCHOR  POINT 
TYPE  10 

10  F0RMAT(/,'  X.Y  COORDINATES  OF  ANCHOR  POINT?  '.$) 

ACCEPT  *.  XA,  YA 

COMPUTE  ANGLES 

DO  100  I  =  1.  471 

IF  Yd)  =YA  .  COMPUTE  THETA  DIFFERENTLY 

IF  (  Yd)  .EQ.  YA  )  GO  TO  20 

FOUR  POSSIBLE  CASES:  XA  >=  XI,  YA  >  YI 

XA  >=  XI,  YA  <  YI 
XA  <  XI,  YA  <  YI 
XA  <  XI,  YA  >  YI 

COMPUTE  THETA  AND  ADJUST  IT  ACCORDINGLY 

THETA  =  ATAN2  (XA  -  X{I),  YA  -  Y(I))  *  Q 

IF  (  YA  .LT.  YI  )  THETA  =  THETA  +  180. 

IF  (  XA  .LT.  XI  .AND.  YA  .6T.  YI  )  THETA  =  THETA  +  3FC. 
GO  TO  60 


lb 


ooo  o  ooo  o  ooo 


c 

20  IF  (  X(I)  -  XA  )  30,  40.  50 

C 

30  THETA  =  90. 

60  TO  60 

(ANCHOR  POINT  AND  LOCUS  POINT  ARE  IDENTICAL:  SET  THETA  =  0) 

40  THETA  »  0. 

60  TO  100 

50  THETA  =  270. 

PERFORM  FINAL  MAPPIN6  OF  THETA 

60  THETA  =  270.  -  THETA 

IF  (  THETA  .LT.  0.  )  THETA  =  THETA  +  360. 

IF  (  THETA  .6T.  360.  )  THETA  =  THETA  -  360. 

100  AN6LE(I)  =  THETA 

SAVE  RESULTS  IN  FILE 

OPEN  (UNIT=1,  NAM£=' LTHETA.DAT',  TYPE='NEW',  ACCESS= ' SEQUENTIAL ' ) 
00  105  I  =  1.  471 

105  WRITE  (1,110)  1+359,  AN6LE(I).  X(I).  Y(I) 

no  FORMAT  (IX,  15,  3X,  F5.1,  2(3X,  F7.5)) 

CL0SE(UNIT»1) 

STOP 
END 
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LAHBOA  THETA  X  T  LAMBDA  THETA  X  Y  LAMBDA  THETA  X  Y  LAMBDA  THETA 
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I  «^^'▼^<bnb')bnbnbnbnbnbnb')bn'0S)'0'C'0'<'O'0'C^^srv^^^vrvr^.r^>rvrs>^<-c^oc^®®XG®o®c>>>•a>c►Q^0'0>->(^0^ 

I  ^s»^v^\^v^■kfM^*»rs^^fs.ls^<rs^vf^^^>^^s.f^.fs.r^.fMr^.rvfM^^f^.^vfVfs.rs.^^.fSkr^.^>.rv^^fv^v^sf^.f^^^vfKf^.^^^v^vr^^^rv^s.^ 

I 

I 


M)Win^^<«fni-)rtrjrtr)f-ifnMnmwfnrnfnmror)Mfnnfnmr>ror)mr)rtmnmr>r)mfownromr)r)«r>r)roMm 
immnvnviinbntnbninintfibiminininininbnmtntnininintninbnminbnv^ininbnintninbninbninininbninininbnbninininin 
«0'0'0'0'0'C'0'0'C'0'<'0«0*0'0'<0'*C'0*0«0^'0*0^'C'0'C'0'0'0'C'0*0'0*0'0'<'4*C'^'0'C'0'0'C'0'0'0'C'0'<<}'0'C'0'<Q 
CM  CM  CM  W  CM  CM  «  CM  fs!  CM  CM  CM  CM  CM  CM  CM  C<  CM  CJ  CM  C'i  C4  CM  CM  Ci  CM  c-i  C'l  C-t  C**  CM  C-t  C-t  C'i  r4  c>:  CM  CJ  CM  Cl  C»  Ci  f  i  Cl  Cl  C:  f4  Ci  C*  Ci  C-t  Ci  Cl  C: 


fnfo«rjwMnmnrjnr)r)mr)Mnrtnwwfnnmmr>mrtfnrtmrtmmfnMwrtmMfnrtr;mmnf*?mmrorofnr>rnm 


I  ooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
iviv)intnininviininininininininininininininininininb'>«ntniov)irbninm;ninmbninininininintnbninbninininbnintin;n 

I  fOfnir'rnrtWfnfnnMfOrtrtfnrjrnfnr)r)#ofO«f*Jrtrninr>rtr>ininrnr>rnr)rnMfnr)rornrofnrtr3fnfnfnr>r’*r5rir;rnr? 


t 

I  o«^CMn^m'orM®OkO«^cMn^vn**orMOO^o^CMro^bnK>rM®o»o^C'im^m'«N®(so^c^;n^bn*crM®a.o^rMr)^ 
I  o»o>o»o^o^o»(S^o»o^oooooooooo^^«^«^^*^w^^«-(CJCMr-4CMCx*C4r’;CMr-^CMr>mmror}mr?mf*;r)^<‘«r^^ 

I 


n^o«Mio<''«CMCM'^««^®rj^»'<cMrv^ri^«^C'imrs.n*«o*br-4rtinrs-*inO'OMO^o*ori©inr-4  0*'Otn'0®-^«c 
p>»®o®®0^ocM^^o^CM^ov)Oinc*in^rsn^«cr»orvinr)wO^i‘vinr<A«4O®rMin^cM«^oors.in^tn«4OO»®r^rM«c 
ir*^mcM»-ood‘0^<<rjin^^«wcMcs^^oc»ooO‘0^o*^ott©a»aNrMrM^rvrsN««'0'«'0'OM)'omininiriin 
®®«o®Q«fNr\rM^fspMfsr^rsf**rscvrvrMfsfMfMfsrv<*^««^'«'0«0'0<'6'<-<)<'«'C'CM)'«««'0'<^'C'0<M>'C«0<' 
CM  CS  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  «  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CN  CM  CM  CM  Cv*  CM  r%  CM  Cv'  CM  CM  CM  CM  CM  CM  M  CM  r-i  CM  CM  CM  <N  CM  CM  CM  CM  CM  W 


n(M^'OCMn^®incM?*)On«^CNO»o^®r)^r>.<cM<®rvnrNO»oo^®Ninn^ino^®CM«co^®CM)no^<‘inM'(M^«r 

o<»i««ib«4^09>CMr:oorMr*?o»^Oo^o^^®CM'OOr)<o«>CM<r'OChocM^'«®o»^ro^mrs.oc«*4M^tf7'Ooo^o^CMCMm 

^»^w«4^«i«*4<MCM<M(M<MCMCMCMCMCMCNCMCMCMCMCccMCMCMr)r«nnmromrororoni*:f*>fnnmr*)K«r;rorcior‘nr*?mnK<r» 


I  *i««4^w<»«^^CMCICMCMrnCMn000000000000000000000000C»0000000000000000 

i»r«mrAtf?®n®®r«®in®®tfnov?iny)inin®ininininininso®v^inv?v)toin®ininvninv)invnv?intninin®»ininm:n 

iwwf*jwinf*i?*?nfnM^fn?orow 

I  »nrorjmr*5rtr?m>omr>wro?or*?rtrjr>r>Wfnmi«inr?r>r>r)f«ronmfOfnrjr#r«fnr? 


im'OfM®CKO»<CM?*iM’«-ors®o*o-<cNrj^*n^fN®o^o*^cMr>^in'^fv®c*o^CMm^in'orv©&^o^cM»r^an'«c«.oo* 

ip:f0«ri?*;^v^v^M*^^«T^inv:viy^inv:bnirwin'0'<'0'<'0'C'«<'0'Cf>»rMr^r^rs?vrs.rvrMfs.©®®ooxm«o« 

I 

I 


o»nno»r9nno^on<^^«oo*4®«4mo«nrNromor)rMV)CM«4CMV^o^'om®^^^^rv©^®omcfrM©r>fMinm'0 
«o®of^^ct^o-«rM^iN^rvin’M‘^c*o^O'0^to^s)OincMOD*i^or>c«»rMOinMrwro^M)^mrMrMrMwoisM'CMo 
I  ^or>r^OM’OCi'00^®CMrs^<<«4^^rsrMoro<Kin«4rs^orMM‘or«»incMCKr\.^CMo®^«'CMOO'rv^^ro(Mo»>®rv 
o«fM^^®^^i-meM*^^ooo*o^©on*rM'0'Cinmnr*^^for)r>cMCMCM^«H*^^*i^ooooco^^»‘^o^C‘^®®® 
^/▼^^^'^^^^♦▼^▼^^^ro*nrt??nMMror>r?r>r>mr?r)r^mr?r)nmr)mfnr)nrjfo«f^CMCMfMrMCMCMCMrMCMrM 


c►fMO^^'00‘fo»•fV'««^v^cMM)M^^®®'or■CMO'Orvr^ro^^.^MWM)^^'OCM-M)®®^v«rnwCK|^•<^•^o*^•-•-o*oCMfnomo 
^o<o^fMO»^cMrM*HO*^cM^<^o<^f^^Ob">^*'^cMCMo*o-^v»!M©®r^»oinc^cM^in^tncMo^'«o*^fv»<inO‘<Mir*c^o^ 
I  CMOkincM®M'  —  fNr»i^yio<<cMf^<M®cMrMCMrs.»^*no^®cino^cM¥n©^Mr^.ocM^rwo^^inio»s©oc-imin'«i^^o*^r-i 
«i««MCMm^?<‘lnbn'CM3rv®©(^c^oo•4^cMc^'m?<?^^5^mlnunM>'C'Or^.r^r^.®®©©ocK^<^CK(>>‘Ooooooo^*«^ 
ininininy)inininninininirtantr’<i<4-M)*<i*0'0'«^'>o«0'0'^'OM>^^M>'M)'0*C'0'<«0'0'0^'«'4^'cr>.rvrvrsrvrvr«.r«rvrv 


I  mininin(ntninm®ininb*)‘/)inviinininv?in!rinr>b')v)b')inin(nb*>v)inv)ininminbnintntnbninininin(Obninmininminio 

I  ?i«?^rn'^fornfor3r5frrnrifn>^rnrr‘r*i»nrifornror>i*:i-ir?rnf«ror^ro»r“rnr:rnfr‘rnr>fn!-jfnforir?ror*,rir;fnr:r>rnr>r; 


I 

I  o««CMnM‘ln'4f^.©l>•o*«^cm^ln'cr<*•X(^o«4r:ro^bn'0^^©c>'0•^r:m«‘bn<>^s©o^c^cMfn^ln'0^«oc>‘0^r'<m^ 
I  ®o®o©©ooo®^CK^Cho-(Kc>>o^>(KOooooooooo*4«-i-«<M^^«4«>«^«Nr-irjrir^r4r^rMr4rMriror)nnro 
I  mininmintn(ninbnb*:x)mbnunbnbntnminbn*CM)'C'0'C'0’<0'M)-^'0'0'0'<}'0'0'0'C>OM)'0'C'0'0-<C'C'0'>0'<}'0'C'0'CM)'0'C 

I 


18 


LAMBIiA  THETA  X  Y  LAHBDA  THETA  X  Y  LAHHOA  THETA  X  Y  LAMBDA  THETA 


>-  I 


I 

I 


I 

I 

I 


) 

I 


I 


I 


I  in  m  b*)  ,r  V) m  »  b*)  <n  bi  ui  b*)  vi  in  in  u?  v)  <n  b'l  u*)  wi  «!  m  m  v)  v:  m 

C4  C-i  C4  C4  C*  C4  Ci  n  C4  (H  C4  C4  Ci  C-J  C  J  C4  C4  C-l  «  C-4  C4  r4  C4  C4  C4  Ci  C4  t4  C4 


ri’-orOfomnr)r)f^'nrtK»ro»^ror>rorjrj»^rOfnMfonfofnmnfnm 

fv.rvrvfsfs.rwfs/s.fs»rs.rvps.ivrs.fs*rsr».rvrsr^rsrsf^rs.r**r**fsisrsrwrs 


t  OOQOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

1  . . . 

i  ininb*)mb‘^inb*)inin(nin(nb*)inbiinininin{ninv?bitnininmirw^(nin 
I  r>»nr>nf^nnr)fnr)rifnMMmmfnr)fnr)fnfOMfnmrjmmfOfom 


O  ^  f  <  -n  ▼  g*)  'O 

o  o  o  o  o  o  o 

O  O  CD  00  C  C  00 


rs  QD  O  »«  r’4 

o  o  o  ^  ^  ^ 

00  O  GO  CD  CD  CO  GO 


<‘m'OrsCDD>o^r4r)^m'Orva)>o 

^^'-»*^-<»<c>4C^r4r^<Me4CMC«C4r4Pi 
OOOOODCDODOOODDOODQO  O  QJ  CD  CD  OD  CD  O 


19 


Appendix  C  -  BASIC  Program  to  Calculate  Wavelength  Needed  for  Mixing 
with  Anchor  Wavelength  to  Produce  Desired  Pair  of  Chromaticity  Coordinates. 


10  REM  THIS  PROGRAM  CALCULATES  WAVELENGTH  NEEDED  FOR  MIXING,  ITS  X.Y 
15  REM  COORDINATES  ON  THE  SPECTRUM 

20  REM  LOCUS  AND  X,Y  COORDINATES  FOR  THE  DESIRED  CHROMATICITY 
40  REM  ACHROMATIC  POINT  IS  SET  AT  D65 

60  REM  IF  YOU  WANT  TO  CHANGE  IT.  ALTER  THE  U'.V  COORDINATES  IN  LINE  80 
80  LET  X0=.1978  LET\Y0=.4683 
90  LET  0=57.2958 

100  PRINT  “ENTER  X  COORDINATE  OF  REQUIRED  DOMINANT  WAVELENGTH; 'MNPUT  XI 
no  PRINT  “ENTER  Y  COORDINATE  OF  REQUIRED  DOMINANT  WAVELENGTH; “MNPUT  Y1 
120  PRINT  “ENTER  THE  FILE  NAME  YOU  WANT  THIS  DATA  STORED  IN  "MNPUT  IS 
130  REM  CONVERTS  1931  CIE  to  1976  UCS 
140  LET  X=4*Xl/(-2*Xl+12*Yl+3) 

150  LET  Y=9*Yl/(-2*Xl+12*Yl+3) 

160  IF  (X-X0)<>060  TO  200 
170  A=90 

180  IF  (Y-Y0)>060  TO  280 
190  A=270S60  TO  280 
200  L=(X-XO)/(Y-YO) 

210  A=ATN(L)*Q 
220  IF  A  >060  TO  250 
230  IF  X>X060  TO  270 
240  A=90-AVG0  TO  280 
250  IF  (X-X0)>060  TO  275 

260  IF  (Y-YO)>OGO  TO  280 

270  A=270-ANG0  TO  280 
275  LET  A=90-A 

300  REM  X9,  Y9  ARE  X.  Y  COORDINATES  FOR  TEST  FIELD 
310  PRINT  “WHAT  IS  X  COORDINATE  OF  DESIRED  COLORT'MNPUT  X9 

320  PRINT  "WHAT  IS  Y  COORDINATE  OF  DESIRED  C0L0R?"MNPUT  Y9 

380  REM  SETS  MONOCHROMATOR  2  AT  ANCHOR  OF  655NM 
390  LET  X2=. 72827/LET  Y2=. 27173 

410  REM  CALCULATES  LENGTH  OF  R2  FROM  655  NM( MONOCHROMATOR)  TO  TEST  POINT 
420  IF  (X9-X2)<>060  TO  460 
430  A2=90 

440  IF  (Y9-Y2)>0G0  TO  550 
450  A2=270\G0  TO  550 
460  L9=(X9-X2)/(Y9-Y2) 

470  A2=ATN(L9)*Q 

480  REM  A2  IS  THETA  FOR  DESIRED  CHROMATICITY  POINT 

490  IF  A2  >  060  TO  520 

500  IF  X9  >  X260  TO  540 

510  A2=90-A2SG0  TO  550 

520  IF  (X9-X2)  >0G0  TO  545 

530  IF  (Y9-Y2)>0G0  TO  550 
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540  A2=278-AA60  TO  550 
545  LET  A2=90-A2 

550  OPEN  “DXOrLTHETA.DAT"  FOR  INPUT  AS  FILE  VFl 
560  OPEN  1$  FOR  OUTPUT  AS  FILE  VF2 
570  FOR  N=1  TO  451 
580  LET  T2=*VFl((4*N)-2) 

590  IF  A2>T2G0>T0  660 
600  IF  A2=T260  TO  620 
610  NEXT  N 

620  LET  L4=VFl((4*N)-3) 

630  LET  X3=VF1  (4*N)-1) 

640  LET  Y3-VF1(4*N) 

650  GO  TO  850 

660  LET  Tl=VFl((4*N)-6) 

670  REM  FINDS  ONE  THETA  ABOVE  T2  AND  STORES  IT  AS  Tl 
680  LET  L2=VFl((4*N)-3) 

690  REM  L2  IS  WAVELENGTH  ASSOCIATED  WITH  THETA 
700  LET  Ll=VFl((4*N)-7) 

710  REM  LI  IS  WAVELENGTH  ASSOCIATED  WITH  THETA  ABOVE 

720  REM  LINEAR  INTERPOLATION  FOLLOWS  FOR  LAMBDA  AND  X  AND  Y  COORDS 

730  L3=(T2-A2)*(L2-L1)/(T2-T1) 

740  REM  L3  IS  AMOUNT  TO  CHANGE  L2  BY  TO  RECEIVE  REQUIRED  LAMBDA 
750  L4=(L2-L3) 

760  REM  L4  IS  LAMBDA  ASSOCIATED  WITH  WAVELENGTH  NEEDED  FOR  MIXING  WITH  ANCHOR 
770  LET  X8=VF1((4*N)-1) 

780  LET  Y8-VF1C4*N) 

790  LET  X7*VFl((4*N)-5) 

800  LET  Y7=VFl((4*N)-4) 

810  LET  X4*  T1-A2)*(X7-X8)/(T1-T2) 

820  LET  X3=(X7-X4) 

830  LET  Y4=  T1-A2)*{Y7-Y8)/(T1-T2) 

840  LET  Y3=(Y7-Y4) 

850  PRINT  "MONOCHROMATOR  WAVELENGTH  FOR  MIXING  IS"L4 

880  PRINT  "X,Y  FOR  MONOCHROMATOR  ARE  "X3,Y3 

882  PRINT  "X,Y  FOR  DESIRED  CHROMATICITY  ARE"X9,Y9 

885  PRINT  “U',V'  COORDINATES  FOR  DESIRED  CHROMATICITY  ARE"X6,Y6 

890  LET  VF2  1)=X3 

900  LET  VF2(2)=Y3 

910  LET  VF2(3)=L4 

920  LET  VF2(4)=X9 

930  LET  VF2(5)=Y9 

940  END 
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